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We consider three- and four-level atomic lasers that are either incoherently (unidirectionally) or 
coherently (bidirectionally) pumped, the single-mode cavity being resonant with the laser transition. 
The intra-cavity Fano factor and the photo-current spectral density are evaluated on the basis of rate 
equations. According to that approach, fluctuations are caused by jumps in active and detecting 
atoms. The algebra is considerably simpler than the one required by Quantum-Optics treatments. 
Whenever a comparison can be made, the expressions obtained coincide. The conditions under 
which the output light exhibits sub-Poissonian statistics are considered in detail. Analytical results, 
based on linearization, are verified by comparison with Monte Carlo simulations. An essentially 
exhaustive investigation of sub-Poissonian light generation by three- and four-level atoms lasers has 
been performed. Only special forms were reported earlier. 

PACS numbers: 42.55.Ah, 42.50.Ar, 42.55.Px, 42.50.Lc 



I. INTRODUCTION 

Interest in the statistics of light emitted by atomic 
lasers has been recently revived as a result of the re- 
alization of micro- lasers [§. The main purpose of this 
paper is to show that expressions derived from rate equa- 
tions coincide with those previously derived from Quan- 
tum Optics. This is so even when the generated light 
exhibits sub-Poissonian statistics. An introduction to 
that method as it pertains to sub-Poissonian light gener- 
ation can be found in tutorial papers ^ . Let us recall 
that light is called sub-Poissonian when the variance of 
the number of photo-detection events over some large 
time duration is less than the average number of events. 
Equivalently, one may say that the spectral density of the 
photo-current is less than the shot-noise level at small 
Fourier (or baseband) frequencies. Analytical expres- 
sions are obtained from rate equations in a straightfor- 
ward manner. Namely, the expression for the internal 
cavity statistics of many 4-level atoms with a negligible 
spontaneous decay previously given in Eq. (4) by Ritsch 
et al. Q is recovered exactly (see Eq. (20) of the present 
paper). Similarly, 3-levels atoms expressions obtained 
by Khazanov et al. are recovered. But coherently- 
pumped 3-level atoms lasers were apparently not treated 



earlier. If the upper and lower decay times of 4-level 
atoms tend to zero, the laser is equivalent to a 2-levels 
atoms laser with Poissonian pump In that limit , the 
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expressions reported in Arnaud and Estcban |7| in |199C 
are recovered. When the above approximations are not 
applicable, the presently reported expressions appear to 
be new. 

The active medium is a collection of identical atoms, 
with the levels labeled |0), |1), |2), and |3) in increasing 
order of energy (see Fig. |l|). Level separations are sup- 
posed to be large compared with /cbT, where T denotes 
the optical cavity temperature and fee the Boltzmann 
constant, so that thermally-induced transitions are neg- 
ligible. Note that the 3-level atoms schemes in Fig. |^(a) 
and mh) are not special cases of the 4-level scheme in 
Fig. ^(c). They are treated separately in appendices. 
Levels |1) and |2) are resonant with the field of a single- 
mode optical cavity. The active medium is supposed to 
be strongly homogeneously broadened so that atomic po- 
larizations may be adiabatically eliminated, an approx- 
imation applicable also to YAG lasers, CO2 lasers and 
semiconductor lasers. Such lasers may operate in the 
steady-state regime, that is, pulsation and chaos do not 
normally occur. 

The probability per unit time that an electronic tran- 
sition from level |1) to level |2) occurs is taken as equal to 
TO, and the probability of an electronic transition from |2) 
to |1) as m + 1, where to denotes the number of photons 
in the cavity (the qualification "per unit time" is hence- 
forth omitted for the sake of brevity). This amounts to 
selecting a time unit whose typical value depends on the 
gain medium. Spontaneous decay from level |2) to level 
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FIG. 1: Level schemes for atomic lasers, (a) A-type 3-level laser, (b) V-type 3-level laser, (c) 4-level atoms laser. For incoherent 
pumping £ — 0. For coherent pumping £ = 1. 



|1) is allowed with probability 7. This decay may be 
either non-radiative or involve radiation into other elec- 
tromagnetic modes, besides the one of interest. Photons 
are absorbed with probability am, where a denotes a 
constant, the absorbing atoms residing most of the time 
in their ground state. These absorbing atoms model the 
transmission of light through mirrors with subsequent ab- 
sorption by a detector. Provided detection is linear and 
reflectionless, it is immaterial whether absorption occurs 
inside or outside the optical cavity. For simplicity, inter- 
nal absorption is neglected. 

"Incoherent" pumping promotes electrons from level 
|0) to level |3) with probability P. When transitions from 
|0) to 1 3) and from |3) to |0) are both allowed with equal 
probabilities, the pumping process is called "coherent" , 
following an accepted terminology. We find it convenient 
to denote by £P the |3) |0) transition probability, with 
£ = for incoherent pumping and £ = I for coherent 
pumping. Coherent pumping is physically realized by 
submitting the atoms to strong optical fields nearly res- 
onant with the |0) — > |3) transition. This pumping field 
may possibly originate from frequency-filtered thermal 
radiation and be highly multimode. Levels |0) and |3) 
need not be sharp. Instead, they may consist of nar- 
row bands for improved coupling to broad-band pumps. 
One-way incoherent pumping would be appropriate to 
describe laser-diode pumps. In semiconductors, however, 
working levels spread into conduction and valence bands, 
so that the present model would not be adequate. Laser 
diodes have been treated previously on the basis of rate 
equations in Spontaneous decay from level |3) to the 
upper working level |2) occurs with probability pu, and 
spontaneous decay from the lower working level |1) to 
the ground level with probability pd, where "u" stands 
for "up" and "d" for down. The relevant probabilities 
are schematized in Fig. |^. 

One of the best-known laser-noise theory is probably 
that of Scully and Lamb ||^ . Incoherent pumping is mod- 
eled by independent injection of 2-level atoms in the op- 
tical cavity. This model leads to a photo-count statistics 



which is, at best, Poissonian. More recently, Ritsch et al. 
Q, Khazanov et al. Ralph and Savage Q consid- 
ered the situation in which the pumping-levels popula- 
tions may fluctuate. At first, it would seem that this 
may only increase the noise. It turns out, however, that 
population fluctuations are correlated in such a way that 
the output light fluctuations may be sub-Poissonian. It 
is difficult to pin point a simple intuitive explanation. It 
has been observed, however, that when lasers are pumped 
through a cascade of intermediate levels, pumping tends 
to be regular |10, a situation somewhat similar to 
laser-diodes high-resistance driving conditions. 

Other means of generating sub-Poissonian light have 
been c onsid ered. Golubev and Sokolov were the 
first in 1984 to point out that lasers with non-fluctuating 



pumps should emit sub-Poissonian light. This conclusion 
has been verified experimentally by Machida et al. jl^ 
with the help of laser diodes driven by high-impedance 
electrical sources. The Scully and Lamb model has been 
generalized to account for regular atom injection p^ . 
Kolobov et al. [|l6| made the interesting observation that 
the photodetection rate spectral density may be below 
the shot-noise level at non-zero Fourier frequencies in the 
case of Poissonian pumps. However, the photodetection 
rate spectral density remains at the shot-noise level at 
zero frequency. Accordingly, such lasers do not generate 
sub-Poissonian light in the sense defined earlier. It has 
been shown that 3-level lasers with coherent decay to the 
ground state generate sub-Poissonian light, and 

further that Raman lasers may generate sub-Poissonian 
light We will not consider here these more exotic 

configurations. A review is in [po[ . 

Rate equations treat the number of photons in the cav- 
ity as a classical random function of time. The light 
field is quantized as a result of matter quantization and 
conservation of energy, but not directly. Rate equations 
should be distinguished from semi-classical theories in 
which the optical field is driven by atomic dipole expec- 
tation values. The theory employed in this paper rests 
instead on the consideration of transition probabilities, as 
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in the Loudon treatment of optical amplifier noise, 
for example. Every absorption event reacts on the num- 
ber of light quanta in the optical cavity. Semi-classical 
theories are unable to explain sub-Poissonian light statis- 
tics because the light generation process and the light 
detection processes are considered separately. The ex- 
pressions obtained from rate equations are found to co- 
incide with Quantum Optics results when the number of 
atoms is large and transitions other than those relating 
to the atom-cavity interaction are incoherent. The laser 
is treated as a birth-death Markov process (see Sec. ||). 
A Monte Carlo simulation gives the evolution of the num- 
ber m of photons in the cavity from which the Fano fac- 
tor ^ = va,r{m)/ (m) is obtained. On the other hand, 
the instants tk when photons are being absorbed provide 
the spectral density of the photo-current, whose normal- 
ized value J5^ is unity for Poisson processes. The nor- 
malized spectrum is denoted The Fourier angular 
frequency f2 is called for short: "frequency". 

When both the number of atoms and the pumping level 
increase, the computing time becomes prohibitively large 
because of the exponentially growing number of events to 
process. Analytical results are obtained by app lying the 
weak-noise approximation. We first give in Sec. Ill A the 
steady-state photon number. The photo-curre nt spe ctral 
density at zero frequency is evalua ted in Sec. 



the photo-current spectrum in Sec. Ill C . The met hod o f 
evaluation of the intra-cavity Fano factor is in Sec. [II D. 
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II. MONTE CARLO METHOD 

The rate-equation model of a A^-atoms singlemode 
laser straightforwardly leads to a master equation for the 
probability of having m photons stored in the cavity at 
time t [|[ ||, |2|] . Alternatively, the laser evolution is mod- 
eled as a temporally homogeneous birth-death Markov 
process. In the steady-sate regime, rate of change and 
equilibrium probabilities of having m and m+1 photons 
within the cavity are linked via a detailed balancing con- 
dition. This is a favorable condition for a Monte Carlo 
simulation |^ because every laser microstate belongs to 
a Markov chain and will thus occurs proportionally to 
its equilibrium probability when the number of step in- 
creases to CXD. 

For the V-type lasers considered in the present section 
(see Fig. |l|(b)), rates of change Wj are ascribed to the 
different kinds of events as given in Table ||. For exam- 
ple, the probability that an atom jumps from level |1) to 
level |3) during the elementary time interval [t,t + St] is 
W2St, where St is chosen small enough that this proba- 
bility be much less than unity. Because atoms are cou- 
pled to one another only through the field, W2 is propor- 
tional to the number ni of atoms in |1) at time t, and 
thus W2 = Pni , where the constant P is proportional to 
the pump strength. If a jump does occur, rii is reduced 
by 1 while the number 713 of atoms in level |3) is incre- 
mented by 1. If the initial value of 7x3 is TV the event 



TABLE L Elementary events in V-type 3-level lasers 
Fig. Wh)) and corresponding rate of change Wj. 



Event 


Transition 


Rate 




photon absorption 






Wi = 


am 


pump absorption 




|3> 


W2 = 


Pni 


pump emission 


13)^ 


|1> 


Ws = 


£Pn3 


coherent emission 


12)- 


|1> 


W4 = 


(m + 1) 712 


coherent absorption 




|2> 


Ws = 


mni 


spontaneous decay 


|2)- 


|1> 


We = 


7712 


upper decay 


|3)- 


|2> 


W7 = 


napu 



does not occur. Similar observations apply to the other 
jump probabilities. Notice that the coherent emission 
rate W4 is proportional to 771 -I- 1, following the Einstein 
prescription. This ensures that laser emission re-starts if 
extinction occurs. A key feature that distinguishes the 
present formulation from other rate-equation methods is 
that absorption of photons by the detector is included 
in the system description. Because detection is supposed 
to be linear, such events are taken to occur with a rate 
Wi = am, where a expresses detector absorption. 

An efficient algorithm has actually been employed : 
Given that an event of any kind occurred at time rj,, the 
next-event time is 



1 



Tk 



In 



(1) 



where r is a random number uniformly distributed in 
the interval [0,1]. The probability that the event is of 
kind I is equal to Wif^- Wj. The Monte Carlo method 
was implemented on a desk computer. Only the total 
number of atoms in each state needs to be tracked (For 
fermions the numerical procedure is significantly more 
involved [|4|). Within the whole time set {t^}, we are 
mostly interested in the subset {tk} of photon-absorption 
events. It is straightforward in principle to evaluate the 
photo-detection noise. We also record m [tk) to evaluate 
the mean value and variance of 777. 

The intra-cavity Fano factor is represented in Fig. ||. 
Both coherent (Fig. |^(a)) and incoherent (Fig. ^(b)) 
pum ping are considered. The analytical results in 
Sec. Ill D , based on linearization agree well with the 
simulation. Notice that is below unity within some 
pumping range. There is good agreement with previous 
Quantum-Optics results p5[ . 

The normalized spectral density is represented 

in Fig. ^ for incoherent (a) and coherent (b) pumps, and 
two sets of parameter- values. For each Monte Carlo run, 
o5^(ri) is first evaluated from the {tk} list |Q and re- 
fined using a smoother power spectral density estimator 
p6| , p7| . Averaging over runs and concatenating neigh- 
boring frequencies produce the final data together with 
error bars at the 95% confidence level. There is fair agree- 
ment between Monte Carlo simulations and analytical 
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FIG. 2: Intra-cavity Fano factor for V-type lasers as a func- 
tion of the pumping rate P. (a) incoherent pumping, (b) 
coherent pumping. Error bars are the 95 % confidence level 
from a statistical treatment applied to ten Monte Carlo runs, 
each havi ng du ration Tm = 200. Plain lines are analytical 
(see Sec. HID). The parameters are: N — 100, p„ = 632, 
a = 6.32. •, 7 = 0; A, 7 = 6.32; 4,7 = 632. 



FIG. 3: Normalized photo-current spectral density of V- 
type lasers as a function of Fourier frequency Q. Points with 
error bars are from Monte Carlo simulations with 150 runs, 
each of d uration Tm = 100. Plain lines are analytical (see 
Sec. QIC). The parameters are: = 100, Pu = 632, a = 
6.32. (a) incoherent pumping •,7 = and P — 1265; 
7 = 6.32 and P = 20. (b) coherent pumping, •,7 = and 
P = 422; A, 7 = 6.32 and P = 20. 



formulas to be subsequently reported. Both predict sub- 
Poissonian photo-current statistics. Even with one bilhon 
photon-absorption events, Monte Carlo spectra exhibit 
large error bars. An analytical method is to be preferred 
when it exists. On the other hand, Monte Carlo simu- 
lations do not rely on linearization and provide a useful 
check. 



III. FOUR-LEVELS ATOMS LASERS 

A. Steady-state 

Let rij, j — 0,1,2,3, denote the number of atoms in 
state j, with 

no + + 712 + n-3 = N. (2) 

The transition rate from |1) — > |2) (stimulated absorp- 
tion) is set equal to the number m of photons in the cav- 
ity, a rule that defines a time unit. Pumping of an atom 
from level |0) to level |3) occurs with rate P. We allow 
for a transition rate £P that an atom in level |3) decays 
back to level |0). Atoms in state |3} spontaneously decay 



to state 1 2) with rate Pu, while atoms in state |1) spon- 
taneously decays to state |0) with rate pd- These spon- 
taneous decays are labelled "upper" and "lower" decay, 
respectively. Note that when 1/pu = the population of 
level 1 3) vanishes and any dependence on £ must therefore 
drop out. Photons are absorbed with rate am. Sponta- 
neous decay from the upper to the lower working levels 
occurs with rate 7. Altogether, there are eight kinds of 
events that may occur in the course of time. 

Let Jj denotes the net pumping rate, TZ the net stim- 
ulated rate, U and V the upper and lower decay rates, S 
the spontaneous decay rate from the upper to the lower 
working levels, and Q the photon absorption rate. The 
steady-state conditions then read 



where 







J = u = v = n 


+ s, 




(3a) 






Q = n, 






(3b) 


J 


^Pno 


- ePns, 


u = 


Puna, 


(4a) 


n 


= (to 4 




V = 




(4b) 


s 


= 7"2, 




Q = 


am. 


(4c) 
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Equations (||), (||) and (||) provide the steady-state 
atomic populations and photon number m. In par- 
ticular 



m 



where 



(5) 

(6a) 
(6b) 

7-1. (6c) 



For moderate pump powers above threshold, m in- 
creases linearly with P. The intercept with the to = 
axis defines the threshold pump power. Mathematically, 
this amounts to first replacing in previous equations m-t-1 
by TO, and then setting to = 0. The laser is found to os- 
cillate for some pump rate provided 



N 




a 




' 1 


2 I +^ 
H 1 
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Vd Pu 








Pd 



7 < (^ - 1) 



+ 1 i + e 



1 -1 



Pd Pu 

The limiting case is obtained at infinite pump power 
(^ - 1) [pd - 7) 



(7) 



lim TO 

P^oo 



(1 



.AM 



-t-2 



■7- 



(8) 



Figure |^ shows how the steady-state photon number 
evolves as a function of the pumping parameter P for 
various values of the spontaneous decay rate 7. Parame- 
ters are similar to |25| except p„-values chosen here for 
a perfect superimposition of the coherent and incoherent 
plots. Moreover it corresponds to optimum pump-noise 
suppression as derived in Se c Jll At very high pump 
levels TO saturates like Eq. (|P because of ground-state 
population depletion. 



B. Zero-frequency noise 

Our analytical results rest on a weak-noise approxima- 
tion. Populations split into steady-state values and fluc- 
tuations. For example the instantaneous photon num- 
ber TO is written as (to) -|- Ato, where (to) denotes the 
steady-state value. Rates split into steady-state values 
and fluctuations consisting of a deterministic function of 
the population fluctuations and uncorrelated Langevin 
"forces". For example J split into J — {J) and A J. The 
latter is the sum of a deterministic function of the pop- 
ulation fluctuations and a Langevin force j{t) expressing 
the jump process randomness. Thus 



J = J + A J, 
n = R + AR, 
S = S + AS, 



U = U + AU, (9a) 
V = D + AD, (9b) 
Q = Q + AQ. (9c) 



m 




FIG. 4: Steady-state number of photons in 4-level atoms 
lasers as a function of P/pd for different values of the sponta- 
neous decay rate "f. N = 10^, a = 6.32, pd = 632. For inco- 
herent pumping Pu = 316. For coherent pumping pu = 949. 
• , 7 = 0; ■, 7 = 6.32; ♦, 7 = 63.2; x, 7 = 632; A, 7 = 6325. 



where 



A J =PAna - f PAng + j, (10a) 
AR = (m + 1) Ani - mAni 

+ {n2 - ni) Am + r, (10b) 

AS" =7An2 + s, (10c) 

AU ^PuAns + u, (lOd) 

AD =PdAni + d, (lOe) 

AQ =aAm + q. (lOf) 

A first-order variation of the expressions in Eq. (^ has 
been performed. 

Conservation of the rates gives 



AJ = AU ^ AD = AR + AS, (11a) 
AQ = AR. (lib) 

Since the total number N of atoms is constant we have 

Ano + Ani + Ana + An^ = 0. (12) 

Replacing atomic populations a nd ph oton numbers by 
their steady-state values in Sec. [II A , the above set of 
equations can be solved. In particular, AQ is a linear 
combination of the Langevin forces 



z£{j^d.u,q.r,s} 



(13) 



where the Cz are real coefficients that depend on the pa- 
rameters N, P, I, Pu, Pd, 7 and a. The detailed expres- 
sions, too lengthy to be given here in their general form, 
are conveniently handled using symbolic calculations. 

The normalized zero-frequency photo-current spectral 
density is of the form 



am ^ — ^ 



(14) 



6 



where cr^ denotes the spectral density value of the 
Langevin noise source z, equal to average rates. 



ar = [m + l)n2 + mni, ad ^ Pdni, 
CTs = 1^2, (Jq ~ am. 



(15a) 
(15b) 
(15c) 



When these expressions are introduced in Eq. (^^ an an- 
alytical expression of is obtained. Three special cases 
are considered below, (a) 7 = and m large compared 
with unity, (b) TV > a, (c) 7 = and TV > a. 

(a) If spontaneous decay is negligible (7 = 0), Eq. (|lj) 
yields 



^ = 1 



, 2(1 



£), 



Pi 



+ 



(6 - 4,yK) ^ 
(^-l)Pd 
2^(2^-pd) 



(16) 



where ^ and ^ are defined in Eq. (^). The normalized 
spectral density is unity at low and high pumping lev- 
els. For some constant value, 5^ reaches its minimum 
valued 



^ min — 

when 

P _ 

Pu 1 



2.yK{^ ~l) + n + i {A-yV {.yV + 15) 



2 (3 + U) {.yV - 1)' 



p _ yV {i + 2£) - {2 + se) 

Vd " (l + ^)(2^-l) 



(17) 



(18) 



(b) When TV > a, Eq. (|T|) yields 

27 4^ + 27 



=5^= 1 



(■?3^ + 7) 



- 7 Pd 

2^(7-Pd)(p<i-2^) 2(l+i?)^F2(7-p^) 



PdPn 



PdPt 



(19) 

Figure ^ gives the normalized zero-frequency photo- 
current spectral density ^ (P jpd-,^ Ivd) in the form of 
contour plots, selecting Pjpu = 2P/pd for the case of 
incoherent pumping and P/pu = ^P/Pd for the case of 
coherent pumping. The darker the area, the lower is the 
spectral density. Since dark areas are wider in Fig. ||(a) 
than in Fig. |5|(b), incoherent pumping is to be preferred. 
Figure |^ shows that sub-Poissonian light generation by 
optically pumped 4-level atom lasers is robust against 
spontaneous decay and pumping conditions. The opti- 
mum conditions (darkest areas) are defined in Eq. (18). 



But small departures from these conditions do not in- 
crease much the noise. 

(c) If both 7 = and N a, the spectral density 
obtained either by setting ^ = cx) in Eq. (^6|) or 7 = 
in Eq. (p^, reads 



y ^ 1 



2PpdPu {Pd + 2{P + Pi + Pu)) 
{2Pp^+Pd{P + P£ + Pu)f ' 



(20) 



TABLE 11: Minimum value of the zero-frequency photo- 
current spectral density =5^mm and intra-cavity Fano factor 
^ for 3 and 4-levels atoms lasers. The conditions on P, pu 
and pd are given. Spontaneous decay from the upper working 
level is neglected and it is assumed that A'" ^ a. 



Laser 






Conditions 


A-type 3-lever 


1/2 


3/4 


Pd = 2P 


A-type 3-level' 


2/3 


5/6 


Pd = 3P 


V-type 3-lever 


1/2 


3/4 


Pu = |P 


V-type 3-level'' 


5/6 


11/12 


P„ = IP 


4-lever 


1/3 


2/3 


Pn = P,Pd = 2P 


4-level'' 


3/7 


5/7 


Pu = 2P, Pd = |P 



"Incoherent pumping. 
''Coherent pumping. 



an expression that coincides with Eq. (4) of |4|. The 
absolute minimum value and corresponding conditions 
are obtained from Eqs. ( p7| ) and ( p^ ) 



5^ 



i + 2e 

3 + Af 



P 

Pu 



1 



p 



l + 2£ 



1 + t Pd 2{l+£) 
£ = 0, 



(21) 



we have therefore 
d = 2P. For coher- 
= 3/7 when p„ = 2P 



For incoherent pumping, i = 
yniin = 1/3 when p^ = P and 
ent pumping, £ = 1, we have 
and Pd = 1^- 

Table O gives the minimum spectral density values 
achievable with optically pum ped three- and 4-level lasers 
as obtained from Eqs. (pi|), (B6) and (A6). Under the 
very special condition of negligible spontaneous decay 
and TV ^ a, the intra-cavity Fano factor depends linearly 
of the zero-frequency normalized photo-current spectral 
density (H), ^ = 2^ - 1. Such relation does not hold 
in the general situatio n whe re Fano factors follow from 
the formulation in Sec 



HID 



V-type incoherently-pumped lasers were treated ear- 
lier by Khazanov et al. Ralph and Savage ^ ex- 
tended the analysis to incoherently pumped A-type lasers 
and 4-level atom lasers. Ritsch et al. gave a descrip- 
tion of 4-level lasers for the two pumping schemes. These 
previous results are exactly recovered from the present 
rate-equations method. 

Formulas derived in the present paper for coherently- 
pumped 3-level lasers, however, appear to be new. Be- 
cause the present method involves simple algebra the con- 
ditions under which the photo-current spectral density is 
minimum are easily obtained. To our knowledge, all the 
results presented in this paper when 7 does not vanish, 
or N is not much larger than a, are new. 



C. Photo-detection spectrum 

The steady-state relations were given in Eq. (pT|). At 
some Fourier frequency the generalized rate equations 
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(a) (b) 




FIG. 5: Contour plots of the zero-frequency normalized 
photo-current spectral density ,5^ for 4-level atoms lasers as a 
function of and P/pd- The laser parameters are the same 
as in Fig. g namely: TV = lO'^, a = 6.32, pd = 632. (a) inco- 
herent pumping, pu = 316, (b) coherent pumping, pu ~ 949. 
In the white area the light statistics is super-Poissonian. 
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FIG. 6: Contour plots of normalized photo-current spectra 
for incoherently pumped 4-level atoms lasers. The laser 
parameters are the same as in Fig. ^ namely: = 10^, 
a = 6.32, pd — 632, p„ = 316. (a) dependence of ^(H) on 
P/Pd with 7 = 0. (b) dependence of on ^/pd with 

P = 316. 



read 



iriAm 


= AR- 


AQ, 




(22a) 


iflAnQ 


= AD- 


AJ, 




(22b) 


iflAni 


= Ai?-|- 


AS- 


AD, 


(22c) 


iflAn2 


= AC/- 


AR- 


AS, 


(22d) 


iflAn^ 


= AJ- 


AU, 




(22e) 



Similar to Sec. |III B| , equations Eqs. (||), (0), ^) and 
( p^ ) are solved for AQ. The formula for the light spec- 
tral density y is the same as Eq. ( p^ ) except that the 
coefficients Cz are complex and frequency dependent 



■ E 

z£{j^d.u.r,s} 



(23) 



where the Langevin "forces" Oz are still given in Eq. (15). 

After rearranging, Eq. ( p3| ) gives the spectral density 
in the form 



,9'iSl) = 1 + 



ao 



where the coefficients Oj and hi are real. 



(24) 



~ 60' 

The form in 

Eq. ( p^ ) ensures that J^(r2) tends to unity (shot-noise 
level) at high frequencies. 

Figure ^ shows that 5^ reaches its minimum value at 
n = 0. When spontaneous decay from the upper work- 
ing level may be neglected, light is always sub-Poissonian 
and the lowest .5^- value occurs when P/pd = 1/2. Spon- 
taneous decay from the upper working level is inconse- 
quential until 7/pd « 310^^. The light statistics ceases 
to be sub-Poissonian when P/pd > 0.3. 



before the number of photons in the cavity. The variance 
of TO equals the integration over frequency of c5^Am(£2), 



the spectral density of Ato. Similarly to Sec. Ill C , 
S^Ami^) obtains by solving for Ato instead of AQ and 
its general form is 



(25) 



where a'^ and &^ are real coefficients. The Fano factor is 
thus 



2n' 



(26) 



The Fano factor for incoherently pumped 4-level atoms 
lasers is represented as a function of pump and sponta- 
neous decay rates in Fig. 0. 

At small pumping rates, and when the spontaneous 
decay rate 7 is large, the photon statistic is essentially 
that of amplified thermal light and the Fano factor ^ = 
(to) -|- 1. The inset in Fig. ^ shows that this simple result 
holds indeed for 7 — 632, excepts for the highest values 
of P, where a slight reduction of occurs. 

When 7 < Pd, the Fano factor exhibits a peak just be- 
low threshold , and decreases to « 1 until the ground 
level gets significantly depleted. Self-quenching |3^ 
is not observed here because atomic polarizations have 
been adiabatically eliminated. Otherwise, our results 
fully agree with those reported by Koganov and Shuker 

i5i. 



IV. CONCLUSION 



D. Fano factor 



The intra-cavity photon statistics is characterized by 
the Fano factor ^ = (Ato^) / (m), where m denotes as 



We have considered optically-pumped 4-level and 3- 
level atom lasers in resonant single-mode cavities. The 
light statistics has been obtained from a simple rate- 
equation approach, using both a Monte Carlo simulation 



8 




with 



FIG. 7: Internal Fano factor of incoherently pumped 4-level 
atoms lasers as a function of P and 7; •, 7 = 0; ■, 7 = 6.32; 
♦ , 7 = 63.2; X, 7 = 632; A, 7 = 6325. The inset is for 
7 = 632. The laser parameters are the same as in Fig. Q 
namely: iV = 10^ a = 6.32, pd = 632. 



and an analytical method based on linearization. The 
emitted light may be sub-Poissonian, as was previously 
observed by many authors. Whenever a comparison can 
be made, exact agreement is reached with the previ- 
ous Quantum Optics results. In the case of coherently- 
pumped 3-level atoms lasers our results appear to be new. 
Note that 3-level atoms lasers are not special cases of the 
4-level scheme. When the assumptions of negligible spon- 
taneous decay and large atom numbers are not made, the 
results presented in this paper for the internal and exter- 
nal field statistics appear to be new. 

For practical reasons, Monte Carlo simulations were re- 
stricted to iV « 1000 atoms. Because the analytical for- 
mulas, obtained through the use of symbolic calculus, are 
lengthy they were not written down in the paper. How- 
ever, they were employed to determine the conditions 
under which the spectral density of the photo-current 
reaches its minimum value. For example, we found that 
when spontaneous decay from the upper working level 
may be neglected, 3-level atoms lasers may deliver light 
with fluctuations at half the shot-noise level. 4-level 
atoms lasers may deliver light with fluctuations at one 
third of the shot-noise level. The photo-current noise de- 
creases further and tends to zero, under ideal conditions, 
when the number of levels becomes large Oj . 



APPENDIX A: A-TYPE LASERS 



The number m of photons in the steady-state is still 



a 



1 



2 + £' 
Vd 



(Ala) 
(Alb) 



^ - 1 - ^ + — ((1 + £) (1 + 7) - '^1) 
Vd 



-1-7. 
Lasing may occur if 



7 < 



.yK - 1 - ^ 



-Vd- 



(Ale) 



(A2) 



^-1-1 

(a) When spontaneous decay may be neglected, 7 — 
0, the zero- frequency normalized photo-current spectral 
density reads 



=y = 1 



2 (1 + f + tJV) 4 (2 -h £) 



Pi 



2[Z + 2l- 2,yV) ,^ 



(A3) 



Note that when ^ — Q this expression may be obtained 
by setting l/p„ = in the 4-level expression. This is not, 
however, a valid procedure in general, 
(b) When iV > a 



27 2(2^ + 7) 



Vd-1 Vd 
{2{2 + £)^ + {4 + £)j) A{2 + i)^'^j 



vl 



v\ 



(c) When 7 = and > a 



.9' = \- 



(A4) 



(A5) 



iPi2 + e)+pdY 

The minimum value of .y and corresponding value of P 
are 

l+i 



9^ 



2 + t 



Vd 
2 + t 



(A6) 



APPENDIX B: V-TYPE LASERS 



The number m of photons in the steady-state is still 



with 



^ = 



N 
2^' 
1 l-f 2£ 



2p« 



— [(2^-l)(7 + 2p„)-l] 



1 7 



^7. 



(Bla) 
(Bib) 

(Blc) 
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Easing may occur if 



(c) When 7 = and iV » a 



7 < 



2^-1 
l + £ + 2£^ 



Pu- 



(B2) 



(a) When spontaneous decay may be neglected, 7 = 
0, the zero-frequency normahzed photo-current spectral 
density reads 



2^ + 1 
(2=yK - 1)^ 



+ 



(1 + 2e) , 

2pS 



(4^ - 1) ^ 
2(2^-l)p„- 



(B3) 



(b) When N :$> a 
27 



^= 1 

^ - 7 Pu 

(1 + 2^) 



^ ^ (1 + 21) - 7) 



+ 



2pS 

2^2^ 



(^_^)(^^ + 2(^-7)p„)' 

(B4) 



^ = 1- 



(P (1 + 2^) + 2p„)' 



(B5) 



The minimum value of ^ and corresponding value of P 
are 



1+41 
2 + AV 



2pu 
l + 2£' 



(B6) 
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